Benchmarked molecular docking integrated molecular dynamics stability analysis for prediction of SARS-CoV-2 papain-like protease inhibition by olive secoiridoids.

Objectives: We performed a virtual screening of olive secoiridoids of the OliveNetTM library to predict SARS-CoV-2 PLpro inhibition. Benchmarked molecular docking protocol that evaluated the performance of two docking programs was applied to execute virtual screening. Molecular dynamics stability analysis of the top-ranked olive secoiridoid docked to PLpro was also carried out. Methods: Benchmarking virtual screening used two freely available docking programs, AutoDock Vina 1.1.2. and AutoDock 4.2.1. for molecular docking of olive secoiridoids to a single PLpro structure. Screening also included benchmark structures of known active and decoy molecules from the DEKOIS 2.0 library. Based on the predicted binding energies, the docking programs ranked the screened molecules. We applied the usual performance evaluation metrices to evaluate the docking programs using the predicted ranks. Molecular dynamics of the top-ranked olive secoiridoid bound to PLpro and computation of MM-GBSA energy using three iterations during the last 50 ps of the analysis of the dynamics in Desmond supported the stability prediction. Results and discussions: Predictiveness curves suggested that AutoDock Vina has a better predictive ability than AutoDock, although there was a moderate correlation between the active molecules rankings (Kendall's correlation of rank (τ) = 0.581). Interestingly, two same molecules, Demethyloleuropein aglycone, and Oleuroside enriched the top 1 % ranked olive secoiridoids predicted by both programs. Demethyloleuropein aglycone bound to PLpro obtained by docking in AutoDock Vina when analyzed for stability by molecular dynamics simulation for 50 ns displayed an RMSD, RMSF<2 Å, and MM-GBSA energy of -94.54 ± 6.05 kcal/mol indicating good stability. Molecular dynamics also revealed the interactions of Demethyloleuropein aglycone with binding sites 2 and 3 of PLpro, suggesting a potent inhibition. In addition, for 98 % of the simulation time, two phenolic hydroxy groups of Demethyloleuropein aglycone maintained two hydrogen bonds with Asp302 of PLpro, specifying the significance of the groups in receptor binding. Conclusion: AutoDock Vina retrieved the active molecules accurately and predicted Demethyloleuropein aglycone as the best inhibitor of PLpro. The Arabian diet consisting of olive products rich in secoiridoids benefits from the PLpro inhibition property and reduces the risk of viral infection.

Molecular Docking and Molecular Dynamics Aided Virtual Search of OliveNet™ Directory for Secoiridoids to Combat SARS-CoV-2 Infection and Associated Hyperinflammatory Responses.

Molecular docking and molecular dynamics aided virtual search of OliveNet™ directory identified potential secoiridoids that combat SARS-CoV-2 entry, replication, and associated hyperinflammatory responses. OliveNet™ is an active directory of phytochemicals obtained from different parts of the olive tree, Olea europaea (Oleaceae). Olive oil, olive fruits containing phenolics, known for their health benefits, are indispensable in the Mediterranean and Arabian diets. Secoiridoids is the largest group of olive phenols and is exclusive to the olive fruits. Functional food like olive fruits could help prevent and alleviate viral disease at an affordable cost. A systematized virtual search of 932 conformers of 78 secoiridoids utilizing Autodock Vina, followed by precision docking using Idock and Smina indicated that Nüzhenide oleoside (NZO), Oleuropein dimer (OED), and Dihydro oleuropein (DHO) blocked the SARS-CoV-2 spike (S) protein-ACE-2 interface; Demethyloleuropein (DMO), Neo-nüzhenide (NNZ), and Nüzhenide (NZE) blocked the SARS-CoV-2 main protease (Mpro). Molecular dynamics (MD) simulation of the NZO-S-protein-ACE-2 complex by Desmond revealed stability during 50 ns. RMSD of the NZO-S-protein-ACE-2 complex converged at 2.1 Å after 20 ns. During MD, the interaction fractions confirmed multiple interactions of NZO with Lys417, a crucial residue for inhibition of S protein. MD of DMO-Mpro complex proved its stability as the RMSD converged at 1.6 Å. Analysis of interactions during MD confirmed the interaction of Cys145 of Mpro with DMO and, thus, its inhibition. The docking predicted IC50 of NZO and DMO was 11.58 and 6.44 µM, respectively. Molecular docking and dynamics of inhibition of the S protein and Mpro by NZO and DMO correlated well. Docking of the six-hit secoiridoids to IL1R, IL6R, and TNFR1, the receptors of inflammatory cytokines IL1ß, IL6, and TNFα, revealed the anti-inflammatory potential except for DHO. Due to intricate structures, the secoiridoids violated Lipinski's rule of five. However, the drug scores of secoiridoids supported their use as drugs. The ADMET predictions implied that the secoiridoids are non-toxic and pose low oral absorption. Secoiridoids need further optimization and are a suitable lead for the discovery of anti-SARS-CoV-2 therapeutics. For the moment, olive secoiridoids presents an accessible mode of prevention and therapy of SARS-CoV-2 infection.